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Abstract 

Intermetallic  compounds  are  under  investigation  as  possible  anode  materials  in  lithium  batteries.  The  compounds  on  which  the  Argonne 
National  Laboratory  group  has  primarily  focused  were  selected  on  the  basis  of  structural  compatibility  between  host  and  product  (lithiated) 
phases.  In  this  presentation,  we  discuss  the  additional  insights  into  the  behavior  of  intermetallic-compound  anodes  that  are  available  from 
materials  modeling. 
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1.  Introduction 

The  presentation  by  Thackeray  at  this  meeting  high¬ 
lighted  measurements  by  the  Argonne  National  Laboratory 
group  on  three  intermetallic  compounds,  Cu6Sn5  [1],  InSb 
[2]  and  Cu2Sb  [3].  The  selection  of  these  compounds  was 
motivated  primarily  by  the  structural  compatibility  of  the 
original  compound  with  the  product  lithiated  phase.  Thus, 
the  sublattice  of  the  more  active  component  of  the  inter¬ 
metallic  compound  (Sn  or  Sb  in  the  cited  compounds)  either 
remains  invariant  or  suffers  at  most  a  moderate  distortion  in 
the  product  phase.  The  concept  of  structural  compatibility 
has  led  to  the  identification  of  a  promising  class  of  materials, 
with  structures  related  to  zinc  blende,  that  have  the  potential 
for  superior  electrochemical  performance.  In  this  presenta¬ 
tion,  we  discuss  the  additional  insights  into  the  behavior  of 
intermetallic  compounds  as  anode  materials  that  can  be 
developed  by  more  detailed  theory  and  modeling. 

In  principle,  modeling  can  address  a  wide  variety  of  issues 
related  to  metallic  electrodes.  The  emphasis  in  our  own  work 
is  on  atomic  structure  and  structural  transformations.  In 
addition  to  the  atomic  scale  phenomena  that  we  focus  on, 
important  processes  occur  at  a  wide  range  of  length  and  time 
scales.  This  multi-scale  aspect  of  the  electrode  kinetics 
makes  their  modeling  particularly  challenging.  We  will 
attempt  to  present  a  viewpoint  that  reflects  a  broad  perspec¬ 
tive,  so  that  the  reader  may  better  understand  what  the 
fundamental  obstacles  are. 
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2.  Structural  compatibility 

If  we  demand  strict  structural  compatibility,  i.e.  that  the 
sublattice  of  the  more  active  element  be  invariant,  and  that  in 
addition  the  sublattice  type  be  cubic-close  packed,  then  the 
number  of  eligible  compounds  is  severely  restricted.  A 
search  of  compound  tables  [4]  reveals  that  only  two  families 
of  III-V  compounds  satisfy  these  criteria.  The  lattice  con¬ 
stants  of  these  compounds,  the  IIIB-antimonide  and  IIIA- 
and  rare-earth  bismuthides,  are  plotted  in  Fig.  1  versus  the 
atomic  number  of  the  cation.  It  is  reasonable  to  suppose  that 
volume  invariance  as  well  as  structural  invariance  is  desir¬ 
able,  and  therefore  also  plotted  in  Fig.  1  (right  ordinate 
scale)  as  horizontal  lines  are  the  lattice  constants  of  the 
lithiated  product  phases  Li3Sb  and  Li3Bi.  The  two  com¬ 
pounds  InSb  and  LaBi  have  the  smallest  lattice  constant 
mismatch,  both  less  than  about  2%. 

InSb  has  been  studied  extensively  by  the  Argonne  group, 
and  by  others,  but  electrochemical  behavior  of  LaBi  has  not 
yet  been  studied.  From  the  standpoint  of  structural  compat¬ 
ibility,  InSb  is  the  most  favorable  system.  As  Thackeray 
pointed  out,  however,  the  performance  of  Cu2Sb  is  superior 
to  InSb,  structural  compatibility  notwithstanding.  Although 
structural  compatibility  appears  important,  it  is,  of  course, 
not  the  only  factor  that  influences  performance.  Some  of  the 
others  will  be  discussed  in  the  following  sections. 

3.  Compounds  under  investigation 

Intermetallic  compounds  have  been  investigated  by  sev¬ 
eral  workers.  Among  the  recently  studied  systems  are  the 
antimonides  Zn4Sb3  [5],  InSb  [2],  SnSb  [6,7],  TiSb2  [8], 
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Fig.  1.  Lattice  constant  data  [4]  for  compounds  of  trivalent  (group  IIIA, 
TUB  and  rare  earth,  RE)  elements  with  Sb  or  Bi.  The  horizontal  lines 
represent  the  lattice  constants  of  Li3Sb  and  Li3Bi.  The  compounds  InSb 
and  LaBi  give  the  smallest  misfit  (<2%)  relative  to  the  displacement 
product  phase  Li3Sb  or  Li3Bi. 

VSb2  [8],  CrSb2  [9],  MnSb  [10],  CoSb3  [11],  Cu2Sb  [3],  and 
CoFe3Sb12  [12],  the  stannides  MnSn2  [6],  Mn!  77Sn  [13], 
Mn3Sn  [13],  Mn3SnC  [13,14],  FeSn  [15],  FeSn2  [8,16], 
CoSn2  [8],  Ni3Sn2  [17],  Cu6Sn5  [1,18,19]  and  CaSn  [20], 
and  the  silicides  CrSi  [21],  NiSi  [22,23],  FeSi  [23]  and 
Mg2Si  [24,25].  There  is  also  some  work  on  chalcogenides 
[12],  particularly  in  earlier  literature  (the  listed  citations 
given  here  are  only  representative,  and  not  exhaustive).  Most 
of  these  compounds  do  not  satisfy  the  criteria  of  product- 
host  structural  compatibility,  and  perhaps  as  a  consequence, 
the  majority  of  them  exhibit  negligible  reversibility  (at  room 
temperature)  of  the  reactions  with  Li  discussed  later.  Thus, 
typically,  the  less  active  component  of  the  compound  is 
extruded  during  the  first  discharge  cycle,  and  the  original 
compound  is  not  restored  upon  cycling.  In  the  later  discus¬ 
sion,  emphasis  will  be  given  to  the  compounds  that  exhibit 
substantial  reversibility. 

4.  Solid  state  reactions  of  lithium 

Three  general  types  of  reactions  of  Li  with  intermetallic 
compounds  are  possible.  The  reaction 

xLi  +  ABy  ->  LixABy  (I) 

is  referred  to  as  a  reconstitution  or  addition  reaction.  In  general, 
the  ternary  product  LixABy  and  the  host  ABy  have  markedly 
different  crystal  structures.  A  second  type  of  reaction, 

xLi  +  ABy  — >  LixBy  +  A,  (II) 

in  which  one  of  the  compound  components  is  extruded,  is 
referred  to  as  a  displacement  reaction.  A  third  type  of  reaction, 

xLi  -I-  ABy  — *  LixAi_zBy  -f~  A  - ,  (III) 


which  we  refer  to  as  mixed ,  involves  both  addition  and 
displacement,  i.e.  only  a  fraction  z  of  component  A  is 
displaced,  and  a  ternary  product  LixAi  _  zBy  is  produced 
as  well.  We  have  adopted  the  convention  in  I— III  that  A  is  the 
less  active  and  B  is  the  more  active  element. 

Displacement  (or  mixed)  reactions  do  not  occur  in  pure 
metals,  but  only  in  compounds.  They  have  previously  been 
studied  primarily  in  oxides  at  elevated  temperatures  [26,27], 
and  little  is  known  about  the  systematics  of  lithium  induced 
displacement  reactions  in  intermetallic  compounds  at  room 
temperature. 

5.  Reaction  sequences 

The  reactions  I— III  represent  the  possible  first  reactions  of 
a  binary  intermetallic  compound.  Typically,  as  the  insertion 
proceeds,  more  than  one  reaction  occurs.  Thus,  e.g.  an 
addition  or  mixed  reaction  may  be  followed  by  a  displace¬ 
ment  reaction  with  respect  to  the  ternary  system,  and  this 
may  be  followed  by  further  addition  reactions  (addition  to 
either  the  host  or  the  phase  of  the  extruded  component,  A). 
We  note  that  for  some  compounds,  the  relevant  sequence  of 
reactions  may  be  ambiguous.  In  this  case,  the  first  principles 
calculations  discussed  below  may  be  helpful  in  elucidating 
the  true  reaction  sequence.  The  calculations  indicate  which 
sequences  are  feasible  energetically,  because  each  reaction 
in  the  sequence  must  have  a  negative  enthalpy. 

6.  Electrochemical  potentials 

In  this  section  we  discuss  the  electrochemical  potential  for 
an  intermetallic  compound  during  the  first  discharge  cycle, 
versus  a  lithium  metal  anode  [28].  To  simplify  the  discus¬ 
sion,  only  the  first  reaction  will  be  considered,  although  we 
have  pointed  out  that  some  compounds  undergo  a  sequence 
of  reactions.  In  Fig.  2,  are  plotted  schematic  potential  curves 
as  a  function  of  lithium  insertion.  The  curve  labeled  ’expt’  is 


Fig.  2.  Schematic  voltage  curve  for  the  first  discharge  of  an  intermetallic- 
compound  electrode  versus  lithium  metal.  The  plateau  region  of  the 
theoretical  curve  corresponds  to  a  two-phase  region,  which  starts  when  the 
solid  solubility  of  Li  in  the  original  compound  is  exceeded.  Since  this 
solubility  is  essentially  negligible,  a  vertical  line  is  drawn  at  zero  insertion. 
Experiment  differs  from  this  behavior  in  that  an  extended  initial  transient  is 
often  observed,  which  we  attribute  to  extrinsic  processes. 
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an  idealization  and  does  not  represent  actual  data  for  a 
particular  system.  The  steeply  sloping  part  of  the  curve  at 
small  insertion  represents  an  initial  transient,  which  we 
attribute  to  extrinsic  processes,  such  as  reaction  with  the 
surface  oxide.  This  transient  is  followed  by  a  nearly  flat 
plateau  region  that  coincides  with  a  reaction  of  type  I,  II  or 
III.  When  the  reaction  is  essentially  complete,  the  electro¬ 
chemical  potential  decreases  sharply.  If  a  sequence  of  reac¬ 
tions  occurs,  then  additional  plateau  regions  (not  shown  in 
the  figure)  may  be  observed. 

The  curve  labeled  ‘theory’  differs  from  experiment  in  that 
the  extended  initial  transient  is  absent.  A  very  narrow 
intrinsic  transient  is  expected  to  occur  as  a  result  of  con¬ 
figurational  entropy,  until  the  solubility  of  Li  in  the  com¬ 
pound  is  reached.  The  equilibrium  solid  solubility  of  Li 
in  most  intermetallic  compounds  at  room  temperature  is 
typically  well  below  0.01.  The  horizontal  line  in  the  theory 
curve  corresponds  to  the  flat  part  of  the  ‘expt’  curve,  but  is 
typically  slightly  above  it  as  shown  later. 


7.  Electrochemical  potential  prediction 


In  a  cell  that  consists  of  a  metallic  Li  anode  and  a  metallic 
cathode,  the  electrochemical  potential  is  defined  as  follows: 


V  =  I  ru  (product)  -  /iLi(Li) 


where  /iLi  is  the  chemical  potential  of  Li  either  in  the  product 
phase  or  in  metallic  Li.  Numerical  calculations  of  the 
electrochemical  potential  for  several  metals  and  compounds 
were  performed  with  the  first  principles  local-density-func¬ 
tional  theory  VASP  code  [29,30],  using  the  generalized 
gradient  approximation  (GGA)  correction.  The  calculations 
correspond  to  zero  temperature,  however  free  energies  at 
room  temperature  differ  only  slightly  from  those  at  lower 
temperatures.  Calculated  electrochemical  potentials  are 
plotted  in  Fig.  3  on  the  abscissa  against  experiment  on 
the  ordinate.  The  calculations  systematically  overestima- 
tethe  voltages,  in  part  because  the  product  phase  is  nucleated 
under  stress  in  the  matrix  of  the  original  compound  or  metal, 
whereas  the  calculations  treat  the  product  as  stress  free. 


8.  Lithium  insertion  kinetics 


Cell  Voltages  (eV) 


Fig.  3.  Experimental  vs.  predicted  cell  voltages  for  Li-induced  reactions 
with  several  intermetallic  compounds  and  pure  metals.  Theoretical  values 
were  obtained  with  the  VASP  code  [29,30],  using  the  GGA.  The 
experimental  atomic  structure,  where  known,  was  employed  as  the  starting 
configuration,  and  lattice  parameter  and  internal  coordinates  were 
optimized.  In  the  case  of  pure  Mn,  an  approximate  antiferromagnetic 
structure  for  the  58-atom  unit-cell  was  assumed.  In  the  case  of  Al2Cu,  no 
electrochemical  reaction  was  observed  experimentally  [8].  Calculations 
were  performed  for  the  ternary  LiAl2Cu,  based  on  a  hypothetical 
hexagonal  structure.  Similar  voltages  were  predicted  for  addition  and 
displacement  reactions. 


and  choosing  a  time  scale  of  5  h,  we  find  that  diffusion  will 
proceed  at  room  temperature  only  if  the  activation  energy  Ea 
is  less  than  about  0.5  eV.  In  relatively  few  metals  have  Li 
diffusion  activation  energies  been  measured.  In  one  system 
of  interest,  InSb,  it  has  been  measured  to  be  0.3  eV,  and  bulk 
diffusion  at  room  temperature  therefore  readily  occurs.  InSb 
is  a  relatively  open  structure,  and  therefore  favorable  for 
interstitial  diffusion.  Whether  bulk  diffusion  at  room  tem¬ 
perature  can  occur  readily  in  other  compounds  is  not  known. 
The  fact  that  the  reactions  with  a  large  number  of  com¬ 
pounds  has  been  observed  indicates  that  the  activation 
energies  either  for  bulk  or  short-circuit  (grain-boundary, 
dislocation,  etc.)  diffusion  is  within  the  range  of  0.5  eV  or 
less. 


The  reactions  described  earlier  must  be  preceded  by  the 
diffusion  of  lithium  into  the  metal  or  compound.  Schema¬ 
tically,  we  can  write  the  diffusion  coefficient  of  Li  inter¬ 
stitials  as  follows: 

D~™2ap(w)J'  (2> 

where  v  is  an  effective  attempt  frequency,  a  is  the  lattice 
constant,  and  kb  is  Boltzmann’s  constant.  For  a  fine-grained 
ball-milled  specimen  with  characteristic  length  /  ~  20  m, 


9.  Volume  expansion 

It  is  well  known  that  one  of  the  problems  with  metallic 
electrodes  is  the  expansion  that  accompanies  the  reactions 
with  Li.  This  expansion  often  leads  to  fracture,  which  may 
result  in  electrical  isolation  of  part  of  the  electrode.  An 
example  is  Si,  which  can  absorb  up  to  4.4  Li  atoms, 

4.4  Li  +  Si  Li4.4Si,  (4) 

and  expands  by  a  factor  of  3.70  (volume  of  product  divided 
by  sum  of  volumes  of  reactants)  in  the  process.  Some  other 
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reactions  are  as  follows: 


3  Li  +  InSb  — >  Li3Sb  +  In 

4  Li  +  InSb  -»  Li3Sb  +  Liln 

2  Li  +  Cu2Sb  — >  L^CuSb  +  Cu 

3  Li  +  Cu2Sb  ->  Li3Sb  +  2  Cu 

10  Li  +  Cu6Sn5  — >  5  L^CuSn  +  Cu 
22  Li  +  Cu6Sn5  — »  5  LLuSn  +  6  Cu 


volume  expansion  =  1 .42 
volume  expansion  =  1.62 
volume  expansion  =  1.51 
volume  expansion  =  1.93 
volume  expansion  =  1.64 
volume  expansion  =  2.84 


These  examples  show  that  there  is  a  wide  variation  in  the 
fractional  volume  expansions  of  different  metals  (or  semi¬ 
conductors),  when  the  reaction  is  carried  to  completion.  On 
the  other  hand,  the  variation  in  the  expansion  per  Li  atom  is 
much  less.  Fig.  4  shows  the  volume  expansion  per  Li  atom 
for  Li-induced  reactions  with  several  intermetallic  com¬ 
pounds  and  elemental  materials.  The  systems  are  plotted 
on  the  abscissa  scale  in  order  of  increasing  volume  expan¬ 
sion.  The  smallest  volume  expansion  is  that  for  LiC6, 
intercalated  graphite,  and  the  largest  is  for  body-centered 
cubic  Li  metal.  Apart  from  InSb,  whose  volume  expansion  is 
similar  to  that  of  intercalated  graphite,  the  other  materials 
plotted  show  a  variation  of  only  about  15%  in  the  volume 
expansion  induced  by  reaction  with  Li.  Not  surprisingly, 
volume  expansion  alone  is  not  a  good  predictor  of  a  materials 
performance  as  a  Li  battery  electrode.  As  discussed  later,  the 
performance  of  an  electrode  is  determined  by  the  complex 
interplay  of  several  processes  that  accompany  Li  insertion. 

The  volume  expansions  presented  earlier  are  based  on  the 
experimental  data.  Calculations  based  on  first  principles 
would  give  very  similar  results.  First  principles  simulation 
can  predict  the  volume  expansion  associated  with  interstitial 
Li  dissolved  in  an  intermetallic  compound,  which  is  difficult 


to  measure  experimentally.  We  present  the  results  of  such 
calculations  here  only  for  InSb  and  Cu6Sn5. 

10.  Kinetics  of  electrode  evolution 

In  the  preceding  section,  we  noted  that  Li  insertion  into  a 
metallic  or  compound  electrode  has  two  immediate  con¬ 
sequences:  (i)  the  equilibrium  solubility  limit  of  Li  is 
exceeded  at  a  low  insertion  level,  and  (ii)  the  Li  induces 
volume  expansion  in  the  material.  Both  the  supersaturation 
of  Li  and  the  induced  volume  expansion  initiate  the  devel¬ 
opment  of  subsequent  materials  processes  that  influence 
intermetallic-compound  electrode  structural  evolution  and 
battery  performance  during  insertion/extraction. 

The  Li  supersaturation  results  in  the  nucleation  and 
growth  of  Li-rich  phases,  in  accordance  with  the  reactions 
I,  II  or  III.  In  the  case  of  reactions  II  and  III,  this  is 
accompanied  by  extrusion  of  the  displaced  component. 
The  volume  expansion  induced  by  Li  results  in  internal 
stress,  plastic  deformation,  rupture  and  spallation.  Thus, 
as  the  insertion  of  Li  proceeds,  atomic  transport  and  che¬ 
mical  reactions  occur,  which  are  influenced  by  internal 
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Fig.  4.  Volume  expansion  per  lithium  atom  for  several  intermetallic  compounds  and  monatomic  metals.  The  host  systems  are  listed  on  the  abscissa  in  order  or 
increasing  specific  volume  of  reaction.  The  volume  of  reaction  corresponds  to  the  difference  between  the  volume  of  the  products  and  that  of  the  original 
compound,  per  Li  atom.  The  volume  of  intercalation  corresponds  to  a  Li  atom  in  an  interstitial  site  in  a  host  matrix.  The  volumes  of  reaction  were  calculated 
based  on  experimental  data,  and  the  volumes  of  reaction  for  InSb  and  Cu6Sn5  were  obtained  from  first  principles  calculations. 
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stresses,  defect  production  and  mechanical  failure.  These 
various  processes  are  all  coupled,  and  cannot  be  treated 
independently.  The  relevant  length  and  time  scales  for  the 
different  processes  range  from  the  microscopic  (e.g.  for 
atomic  diffusion)  to  the  macroscopic  (e.g.  for  fracture  and 
spallation). 

To  model  electrode  performance,  in  principle,  would 
require  the  development  of  a  set  of  kinetic  equations,  which 
ideally  would  incorporate  all  of  the  relevant  physical  (che¬ 
mical,  atomic-transport  and  mechanical)  processes.  The 
solution  to  these  spatio-temporal  equations  for  suitably 
chosen  boundary  conditions  (such  as  external  current  and 
geometry)  would  give  information  about  the  evolution  of  the 
electrode  during  battery  operation.  The  voltage  profile  as  a 
function  of  time  and  the  capacity  would  be  determined,  as 
well  as  information  about  evolution  of  structure,  phase 
distribution,  and  composition.  Discharge  and  charge  cycles 
would  be  treated  by  simulating,  in  turn,  external  currents  in 
opposite  directions.  Capacity  retention  could  be  explored  by 
examining  the  behavior  upon  repeated  cycling. 

Unfortunately,  the  battery-operation  phenomena  that  we 
wish  to  describe  go  beyond  the  current  state  of  the  art  of 
kinetic  modeling.  The  most  significant  technical  obstacles  to 
modeling  the  kinetics  of  the  intermetallic  electrodes  are  (i) 
the  multiple  length  and  time  scales  involved,  mentioned 
earlier,  and  (ii)  the  inhomogeneity  of  the  atomic  transport 
and  structure  evolution  (the  activation  energy  for  bulk 
diffusion  likely  exceeds  the  value  of  0.5  eV,  the  approximate 
upper  limit  for  processes  to  be  active  at  room  temperature). 
Other  complications  are  the  behavior  of  the  solid-electrolyte 
interface  [31],  which  undergoes  continuous  evolution  during 
battery  operation,  and  the  binder  [32]. 

The  existing  methodologies  for  kinetic  modeling  include 
first  principles  theory  [29,30]  and  molecular  dynamics  [33], 
which  are  restricted  primarily  to  microscopic  length  and 
time  scales,  Monte  Carlo  lattice  gas  models  [34],  which 
extend  to  mesoscopic  length  and  time  scales  and  continuum 
and  rate-theory  models  [35,36],  which  extend  to  macro¬ 
scopic  length  and  time  scales.  Continuum  and  rate-theory 
models  are  widely  practiced  in  kinetic  modeling.  For  exam¬ 
ple,  they  are  commonly  used  to  model  diffusion  limited 
growth  of  precipitates  [37].  Electrochemical  applications, 
compared  with  precipitate  growth  kinetics,  have  several 
additional  complications,  including  external  current,  rupture 
and  spallation,  simultaneous  diffusion  of  different  compo¬ 
nents,  diffusion  governed  by  short  circuit  paths  (planar  and 
linear  defects),  the  SEI,  etc.  These  complications  have  thus 
far  discouraged  attempts  to  model  intermetallic  compound 
battery  electrode  kinetics. 

11.  Conclusions 

We  have  argued  earlier  that  the  complexity  of  the  materi¬ 
als  processes  and  transformations  that  occur  during  battery 
operation  make  performance  simulations  for  intermetallic 


compound  electrodes  intractable  at  this  time.  Nevertheless, 
in  spite  of  our  inability  to  model  the  overall  kinetic  behavior 
of  the  electrode,  some  thermodynamic  and  microscopic 
aspects  of  intermetallic  electrodes  can  be  modeled  with 
reasonable  reliability.  These  properties  include  the  reaction 
voltage,  discussed  earlier,  which  first  principles  calculations 
typically  overestimate  slightly,  but  give  the  correct  order  of 
magnitude  for.  Such  energy  calculations,  in  conjunction 
with  in  situ  structure  measurements  such  as  X-ray  diffrac¬ 
tion,  can  also  help  identify  the  actual  sequence  of  reactions 
of  Li  with  an  intermetallic  compound  that  occur  during  the 
insertion  cycle.  Another  property  that  is  accessible  to  theory 
is  the  partial  molar  volume  for  Li  insertion  into  a  metal. 
Theory  can  predict  properties  of  intercalated  (interstitial) 
lithium,  such  as  site  preference,  solubility  and  perhaps 
activation  energy  for  diffusion,  although  such  calculations 
were  not  presented  here.  Therefore,  although  the  overall 
reaction  kinetics  are  beyond  our  present  ability  to  model, 
theoretical  predictions  can  be  helpful  in  the  interpretation  of 
experiments  on  particular  intermetallic  compounds. 
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